Quantum dot with ferromagnetic leads: a densiti-matrix renormalization group study 
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A quantum dot coupled to ferromagnetically polarized one-dimensional leads is studied numer- 
ically using the density-matrix renormalization group method. Several real space properties and 
the local density of states at the dot are computed. It is shown that this local density of states is 
suppressed by the parallel polarization of the leads. In this case we are able to estimate the length 
of the Kondo cloud, and to relate its behavior to that suppression. Another important result of 
our study is that the tunnel magnetoresistance as a function of the quantum dot on-site energy is 
minimum and negative at the symmetric point. 

PACS numbers: 73.63.Kv, 72.15.Qm, 72.25.-b 



The relentless pursuit of smaller, faster, more efficient 
electronic devices has led to an increasing interest in 
mesoscopic and nanoscopic devices behaving as quantum 
dots (QD), using semiconductor technologjii and nan- 
otubes as components^ among other possibilities. Some 
of these devices allovif the exploitation of the electron spin 
in addition to its charge with various potential applica- 
tions in the so-called spintronics.^ In particular, a spin- 
tronic device connecting a QD to ferromagnetic leads has 
already been experimentally studiedi^ 

These nanoscopic devices have also been growing in 
complexity. For example a device consisting of two QDs 
has been developed to show a nonlocal control of one QD 
over the other via a RKKY interactionn^ Another exam- 
ple is the device where an appropriate arrangement of 
electrodes defines a section of a nanotube as the QDf& 
This QD turns out to be then connected to essentially 
one-dimensional leads. In the first example, a real space 
description of the magnetic correlations would be impor- 
tant to understand the internal working of this device. 
In the second example, the effects of finite length of the 
nanotube sections connected to the QD would have to be 
considered. 

Among many techniques used to study these kinds of 
problems, which have been developed in the context of 
the Kondo effect,^ the numerical renormalization group 
(NRG) has provided many important and essentially ex- 
act results. However, this technique cannot tackle the in- 
creasingly complex devices such as the above mentioned 
ones. For this reason, we propose in this paper the appli- 
cation of a real space numerical technique, the density- 
matrix renormalization group (DMRG)*, which has been 
extensively used to study quasi-one-dimensional strongly 
correlated electron systems i2i This technique is ideally 
suited to provide detailed real space information, such 
as spin-spin correlations or electron site occupancies, re- 
quired to understand nanoscopic devices. In addition, it 
works on finite size systems. This is not a disadvantage 
when bulk leads are involved in the device because in 
most cases one could resort to some kind of extrapola- 
tion. On the contrary, it can provide results for finite 
systems such as the example mentioned before. 



In this paper, we will study the simplest device in spin- 
tronics, the spin valve. This device consists of a single- 
level QD attached to two noninteracting leads. On these 
leads, a ferromagnetic polarization is introduced by an 
applied magnetic field, or in other words, by a rigid dis- 
placement of spin-up and spin-down electron bands. The 
central quantities to study are the conductance when the 
polarization of the leads are parallel (Gp) and antiparal- 
lel (Gap)- The main measure of the performance of the 
spin valve is the tunnel magnetoresistance (TMR) defined 
as TMR = [Gp— Gap) /Gap- This device has been stud- 
ied before using different approximationsi2iii*iS. These 
analytical techniques are not exact and in fact they lead 
to some degree of controversy around central issues such 
as the behavior of the conductance when the leads have 
a parallel polarization. Some of these controversial issues 
have been partially settled using NRG>i2*i^ 

As in previous studies on this system, the Hamiltonian 
model for a QD located at the center of the chain is a 
generalization of the Anderson model defined as 

- t J2 {clc.+ia+H.c.)-hn ^ S: 

i=>l,<T i=>l 

- t'J2i(^-laC0a + cl^Ci^ + H.C.) 

a 

-I- e no + C/no,tno,j, (1) 

where conventional notation was used. We adopt t as 
the scale of energy. The magnetic field on the leads 
can be set in two configurations, hj^ = hj^ = h, which 
corresponds to a parallel (P) polarization of the leads, 
and = —hji = h which corresponds to the antiparal- 
lel (AP) one. In the following, h > favors a positive 
polarization. Notice that the leads are described by a 
real-space tight-binding model which is in principle more 
realistic than the ones implicit in the NRG treatment. 

Model equationQJ was studied by the DMRG on 
L — 63, 79 and 95 chains with open boundary condi- 
tions. Most of the results shown below were obtained 
for two sets of parameters: U = 1, t' — 0.4 and U = 4, 
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FIG. 1: (Color online) The LDOS for U = 1, t' = 0.4, and 
e — —0.5, for various values of the polarization p, (a) AP 
magnetization for the L — 79 chain, and (b) P, L = 95. The 
inset shows the length of the "compensation cloud" (see text) 
for the P case, L — 95. 

t' = 0.8. For both sets of parameters the effective Kondo 
coupling at the symmetric point is J = 4i'^/C/ = 0.64, 
although strictly speaking this relationship is only valid 
for U » t' . There are, however, many important dif- 
ferences between both sets of parameters as we will show 
below. 

The main quantity that we have studied is the lo- 
cal density of states (LDOS), p{uj), at the QD. In the 
first place, from this quantity it is possible to evaluate 
the conductance in the linear response regime. In the 
second place, recent advances in scanning tunneling mi- 
croscopy (STM) have made it possible to directly mea- 
sure this quantity. Since in our DMRG calculation the 
ground state vector is measured when the two added sites 
are at the center of the chain (symmetrical configura- 
tion), the QD then is one of these two sites which are 
exactly treated. Then, we adopt the approximation of 
applying the creation and annihilation operators at the 
QD on the ground state vector and then determine p{uj) 
following the well-known continued fraction formalism. 
A more accurate approach would be, after the applica- 
tion of each of those creation and annihilation opera- 
tors, to run additional sweeps for an enlarged density 
matrixii^ In any case, the truncation of the Hilbert space 
is the essential source of error in the DMRG, and to es- 
timate the precision of our approach we have compared 
results for various numbers of retained states A/, from 
M = 500 to 1000. We would also Hke to stress the 
fact that the conductance in linear response is related 
to the LDOS near lu — where the approximation is 
more precisesiSi In addition, we have computed standard 
properties such as the electron occupancy of each site, 
< rii o- > (cr =tj J.) and spin-spin correlations from the 
QD, 'S{j) =< S§ S^>-< S§ X >. 

All the results presented correspond to the half-filled 
system. For the AP case, we work in the = 1/2 sector. 
For the parallel polarization, we work in the S''^ subspace 




FIG. 2: (Color online) The LDOS for U = 4, t' = 0.8, 
and e = —2.0, for various values of the polarization p, (a) 
antiparallel, L = 79 (solid lines), L = 95 (dashed lines); 
(b) parallel configuration, from top to bottom at w = 0: 
L = 79, p = 0.013,0.203,0.321,0.426 (solid lines); L = 95, 
p = 0.011,0.399 (dashed lines). The inset shows the length of 
the "compensation cloud" (see text) for the P case, L — 79 
(open circles), L = 95 (filled circles). 

which minimizes the total energy, and hence it depends 
on h. The polarization of each lead Pa (a = L, R) is 
defined as Pa — (jt-qi — nai)/{na'\ + Uai). Then, the 
polarization p of the P (AP) configuration is {pl +pr)/2 
[{pl ~Pfl)/2]. For small clusters, and specially for the 
U — 1 case, it is difficult to get pl — pr, i.e. the same 
numbers of solitons on each lead. In these cases the AP 
case shows spurious suppression of the LDOS.^« 

In Fig.Qwe show the LDOS at the QD, {/ = 1, i' 0.4, 
at the symmetric point e = — 1//2, for several values of 
the polarization p. In this figure and all the following 
similar ones, we adopted a Lorentzian broadening of the 
peaks of (5 = 0.1 In the first place, for p = 0.013, cor- 
responding to /i = 0, the LDOS is the one expected in 
the Kondo problem. The small splitting of the Kondo 
peak is a finite size effect. For the case of an AP con- 
figuration (Fig. nja)), a small value oi p — 0.088 leads 
to a slight modification of the LDOS, that virtually re- 
mains unmodified by further increasing the polarization 
of the leads, at least in the range examined {p < 0.3). 
The largest value considered, p = 0.3, is obtained for 
h — 2.2 {L — 95). On the other hand, for the P con- 
figuration (Fig. n^b)) it can be seen as a suppression of 
the Kondo peak for the polarization p > 0.2. For small 
values of p the peak at w = remains roughly unmodi- 
fied. As we will show below, this different dependence of 
the LDOS with p for the AP and P arrangements leads 
to the expected behavior of a spin valve. The inset in 
Fig. 1 shows, for the P magnetization, the length of 
the "compensation cloud" i^, a possible measure of the 
Kondo cloud, defined in such a way that the sum of the 
correlations S{j) (— C/2 < j < £,/2,j ^ 0) is equal to 
0.9S'(0)ii The reduction of ^ with p, due to the better 
screening of the minority spin at the QD by the majority 
spins on the leads, is consistent with the suppression of 
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FIG. 3; (Color online) Sf for the L = 79 chain, U ^ 4, t' ^ 
0.8, e = -2.0, for p fa 0.01,0.1,0.2 and 0.3, from bottom to 
top on the left lead; (a) antiparallel, (b) parallel configuration. 
The inset shows 5^ at the QD as a function of the polarization 

the Kondo resonance observed in the LDOS. 

Qualitatively the same behavior can be observed for 
the case of the parameters U — 4, t' = 0.8, at the sym- 
metric point e = —2.0. In this case, the Kondo resonance 
has a smaller weight as compared with the previous set of 
parameters. In Fig.|2a), which shows the LDOS around 
w = 0, it can be seen again that for the AP case the LDOS 
is not much changed as p is increased. In contrast, for 
the P case, the Kondo peak is strongly suppressed in the 
presence of polarization. This suppression is considerably 
stronger than for the previous set of parameters, U = 1, 
t' = 0.4. We have also examined the size dependence 
of these results by computing the LDOS on the L = 95 
chain. Although the position of the peaks, as expected, 
moves to lower frequencies (in absolute value) , the overall 
behavior with p is similar to that for L — 79. As in the 
previous case, the length of the "compensation cloud" 
(shown in the inset in Fig. 2) for the P polarization 
decreases with increasing p in agreement with the sup- 
pression of the Kondo resonance. It is interesting to note 
that, consistently with this behavior, the value of ((S'o)^) 
also decreases from ^ 0.18 for p ^ to 0.10 for p ^ 0.4. 
Notice that in the full Kondo regime {{S§y) ^ 0.25. 

The mean value of is shown in Fig. |3| for J7 = 4, 
t' — 0.8, and e — —2.0, for the AP and P polarizations 
in panels (a) and (b), respectively. The inhomogeneous 
way in which the sites are polarized, due to the open 
boundary conditions of our clusters, is clearly visible. We 
have also observed this behavior for the same parameters 
on the 95-site cluster, suggesting that it is not a mere 
finite size effect. The different behavior of the AP and P 
configurations can be traced to the behavior of the mean 
value of the spin at the QD. The apparently different 
behavior of Sq — Sqj^ as a function of the polarization 
is shown in the inset. For the AP arrangement is 
almost not affected by p, while for the P arrangement 
S'o goes to the minimum possible value by increasing p. 
This behavior of in the P case is well known and 
it has been shown that a magnetic field applied to the 
QD could in fact reduce the suppression of the LDOS 
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FIG. 4: (Color online) The LDOS for the L = 79 chain: U = 
1, t' = 0.4, and e =-0.5 (solid line), -0.25 (dotted line), 0.0 
(dashed line), 0.25 (dotted-dashed line), 0.5 (dotted-dotted- 
dashed line) for (a) p = 0.0, (b) p = 0.30, AP configuration, 
and (c) p = 0.30, parallel configuration. 

discussed ahovem^^ 

We now study the most important issue in this kind of 
systems, i.e., the behavior of the conductance and hence 
of the TMR as a function of the applied voltage. To un- 
derstand this behavior we show in Fig.0|the LDOS on the 
L = 79 chain for the parameters U = 1, t' = 0.4, lead po- 
larization p ~ 0.3, and various values of e. In Fig.^a) we 
show for comparison the unpolarized case. The relevant 
quantity is the LDOS at lu — which is proportional to 
the conductance, G = {2Trt'^ / 1) p{uj = 0), in units of e^/h. 
G as a function of e has the typical form of a dome as it 
can be seen in Fig.jSJa). However, p{oj — 0) for a finite 
cluster is an artifact of the Lorentzian broadening S of the 
peaks adopted in our calculations. Although p{lu = 0) 
could be computed in this waj*i^, we prefer to discuss 
the results of Fig. ^ at a qualitative level. The LDOS 
for p = 0.3 in the AP polarization is shown Fig. 0{b). 
Its to dependence for various values of e is very similar 
to the one for the unpolarized case. On the other hand, 
for the P configuration, the LDOS, depicted in Fig.^Jc), 
shows at a; = a nonmonotonic behavior as a function 
of the gate voltage. This nonmonotonic behavior, when 
translated to the conductance, is crucial to understand 
the e dependence of the TMR. 

Our most relevant results are shown in Fig. |31 To 
compute the conductance we use the Friedel sum rule^, 
Gcr = sin^ nria, where is the occupancy at the QD of 
an electron with spin a. The Friedel sum rule is valid for 
arbitrary polarization in the P case, except for negligi- 
ble higher order correctionsi^. The total conductance is 

FigureEIa) shows the conductance for the unpolarized 
case, AP and P polarizations, for U = 1, t' = 0.4, and 



4 




I I I I I I I I I I I L 

-5 -4 -3 -2 -1 
8 



FIG. 5: (Color online) (a) Conductance as a function of the 
gate voltage for U = 1, t' = 0.4, L = 79, and p = (open 
circles). Results for L = 63, p = 0.30: P configuration (open 
triangles), TMR (filled triangles), and L = 95, p = 0.35 AP 
(crosses), P (plusses), and TMR (stars), (b) Same as (a) for 
U = 4, t' — 0.8. Symbols for L = 63 and L = 95 same as (a). 
Results for L = 79: AP (squares), P (diamonds), and TMR 
(filled circles). The values of the polarization are p — 0.31 
(L = 63), p = 0.32 {L = 79), and p = 0.36 (L = 95). (c) Up- 
spin (circles) and down-spin (squares) contributions to Gp for 
the parameters of (b), L = 79. 

Fig. Efb) for [/ = 4, t' = 0.8. The conductance for the 
AP configuration is virtually indistinguishable from the 
unpolarized one for L — 79. On the other hand, Gp has 
a minimum at the symmetric point, then increases as e 
moves apart from that point, and finally it merges with 
the curve for the unpolarized case, decreasing for large 
values of the gate voltage. Then, in the most important 
region for application purposes, i.e. near the symmetric 



point, the TMR is negative with a minimum precisely 
at this point. This result is in contradiction with the 
one reported in Ref. Q, showing a positive value at the 
symmetric point. This difference may be traced to the 
model we adopted for the leads, which implies loosing the 
particle-hole symmetry for the P case, which is present 
m Ref.Eisa Notice that for the case of Fig.Etb), there 
is a clear trend as the chain size is increased, that is, 
a deepening of Gp around the symmetric point, with 
a subsequent decrease in the TMR. Notice that the lead 
polarization is slightly different in the three clusters. This 
trend is not that clear for the case of Fig.|SJa), since here 
the smaller value of U requires larger lattice sizes. 

Finally, the up- and down-spin contributions to Gp are 
shown in Fig. |5{c) for the couplings of Fig. E^b) . At the 
symmetric point G| = G|. As e is moved away from 
the symmetric point, the conductance of the minority 
spin becomes dominating, reaching perfect conductance, 
at e ~ —U and e ~ 0, and it accounts for the increase of 
the total conductance. 

In summary, we have shown the possibility of apply- 
ing the DMRG technique to the problem of transport 
through quantum dots. For the particular problem of a 
quantum dot coupled to polarized leads we have shown 
that the LDOS at the QD is suppressed by the P polar- 
ization of the leads. In addition, we have shown that the 
TMR as a function of the quantum dot on-site energy 
is minimum and negative at the symmetric point, a re- 
sult which is at variance with those reported in previous 
studies. This is a central issue in these kinds of devices. 
In addition to this, we have shown that this technique al- 
lowed us to estimate the length of the Kondo cloud, and 
to relate its behavior with the suppression of the Kondo 
peak, correlated in turn with the conductance. This em- 
phasizes the fact that real space properties, which are 
accessible by the DMRG, and are very difficult to obtain 
by alternative techniques, are very important to under- 
stand the behavior of these devices. 
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